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Abstract 

We compute 1-loop corrections to Lorentz-signature de Sitter-invariant 2-point functions defined 
by the interacting Euclidean vacuum for massive scalar quantum fields with cubic and quartic 
interactions. Our results apply to all masses for which the free Euclidean de Sitter vacuum is 
well-defined, including values in both the complimentary and the principal series of SO(D, 1). In 
dimensions where the interactions are renormalizeable we provide absolutely convergent integral 
representations of the corrections. These representations suffice to analytically extract the leading 
behavior of the 2-point functions at large separations and may also be used for numerical compu- 
tations. The interacting propagators decay at long distances at least as fast as one would naively 
expect, suggesting that such interacting de Sitter invariant vacuua are well-defined and are well- 
behaved in the IR. In fact, in some cases the interacting propagators decay faster than any free 
propagator with any value of M 2 > 0. 
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I. INTRODUCTION 



While free quantum fields in de Sitter space (dS) have been well understood for some time 
(see [1] for scalar fields), interacting de Sitter quantum field theory continues to be a topic 
of much discussion. In particular, the literature contains numerous suggestions of possible 
quantum field theoretic instabilities (see e.g. [2-5]), many of which have been argued to 
perhaps lead to decay of the effective cosmological constant. Our goal here and in [6] is to 
address the specific class of such concerns associated with infra-red (IR) divergences of the 
naive Lorentz-signature de Sitter Feynman diagrams, or more generally those concerns that 
can be addressed in the context of minimally-coupled scalar fields with M 2 > 0. 

As we will review, IR divergences arise in generic scalar field theories in Lorentz-signature 
perturbation theory about the free Hadamard de Sitter-invariant vacuum. (This vacuum 
is often called the free Euclidean vacuum as it may be defined by analytic continuation 
from Euclidean signature.) While such divergences can be avoided at tree level when the 
fields are sufficiently heavy, they nevertheless arise in loop diagrams. On the other hand, 
since Euclidean de Sitter is just a sphere, it is clear that there are no IR divergences in 
Euclidean signature. Our goal is to demonstrate that that no pathologies arise from ana- 
lytic continuation of interacting Euclidean vacuua to Lorentz signature, where they define 
de Sitter- invariant states. Specifically we show that, at least through 1-loop order, the as- 
sociated Lorentz-signature 2-point functions for massive scalar fields with cubic and quartic 
interactions are finite and decay at large separations at least as fast as one would naively 
expect. This indicates that these Lorentz-signature de Sitter invariant vacuua are both 
well-defined and well-behaved in the IR. In particular, it suggests that these vacuua are 
stable. 

Our results apply to all masses for which the free Euclidean de Sitter vacuum is well- 
defined, i.e. for all M 2 > 0, including values in both the complimentary series and the 
principal series of SO(D, 1). In dimensions where the interactions are renormalizeable, we 
provide absolutely convergent integral representations of the corrections which allow us an- 
alytically extract the leading behavior of the 2-point functions at large timelike separations. 
In addition, the representations are amenable to numerical calculations, demonstrating that 
our methods provide practical tools for calculating Lorentz-signature correlation functions. 
We provide a number of checks on our results, including consistency with known flat-space 
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limits. The complications associated with both higher loops and higher n-point functions 
will be addressed in [6], with similar conclusions. Such results are in qualitative agreement 
with those obtained using stochastic inflation techniques [7] , which are expected to be valid 
in the limit Mi <C 1, where i is the de Sitter length scale. 

We begin by briefly reviewing de Sitter field theory in section II, and by reviewing some 
useful tools for analytic continuation in section III. We then compute perturbative correc- 
tions to propagators in section IV and establish their IR properties, though some details 
are relegated to the appendices. An interesting feature is the fact that, in some cases, the 
corrections enhance the fall-off of the propagator at large times by opening what is effec- 
tively a decay channel, even when the daughter particles are heavier than the field under 
consideration. This corresponds to the fact that particles in de Sitter space can decay to 
heavier particles (see e.g. [8-11]) due to the lack of a globally timelike Killing vector field 
(so that there no conserved notion of energy that is positive definite). When this occurs, 
the fall-off of the corrected propagator can be faster than that of any free field with M 2 > 0. 
Section V then closes with a summary and discussion of general stability issues for de Sitter 
space. 

II. FIELD THEORY IN DE SITTER SPACE 

The following brief review of de Sitter scalar field theory provides an opportunity to fix 
conventions and to discuss the IR divergences of naive de Sitter Feynman diagrams. Consider 
the D-dimensional de Sitter space dSo for which the metric in global coordinates is 



where dVl 2 d is the metric on the unit d — D — 1 dimensional sphere S d . Free scalar fields 
obey the Klein-Gordon equation 



and define representations of the (connected) de Sitter group SOq{D, 1). It is useful to 
define the dimensionless mass parameter a by —a (a + d) :— M 2 i 2 . Throughout most of our 
work, the ambiguity a — > — (a + d) will be a redundancy of our description, and symmetry 
a — > —(a + d) will provide a useful check on our calculations. However, for the moment 



ds 2 = i 2 [-dt 2 + (cosht) 2 dn 2 d ] , 



(1) 



n x( j)(x) = M 2 <p(x), 



(2) 
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choosing the branch 

1/2 

(3) 

the standard de Sitter representations may be classified as follows [12]: 

1. complementary series: — d/2 < a < , 

2. principal series: a = —d/2 + ip, p G R, p > , 

3. discrete series: a = 0, 1, 2, . . . . 

We denote the line Re a = —d/2 on which the principal series fields live by Tp below. 

Fields with M 2 > correspond to representations in the complementary and principal 
series (see figure II). In particular, sufficiently light fields belong to the complementary 
series while heavier fields belong to the principal series. In either case, a de Sitter-invariant 
Green's function A° y (with arguments x,y) may be defined by analytic continuation from 
Euclidean signature (i.e., from S D ). We summarize this construction in section III below, 
but for now we merely state that in the principal and complimentary series the propagator 
A°" contains terms that fall off like e CT '*', e~^ a+d ^ when one argument is held fixed and the 
other is taken to large values of \t\. It is important to note that the fastest such decay occurs 
in the principal series where Re a = Re[— (a + d)] = —d/2. We will not consider massless or 
tachyonic scalars further, as the corresponding free theories do not admit de Sitter-invariant 
Green's functions [1]. 

Let us now briefly review the IR diverges that arise in calculating naive Lorentz-signature 
Feynman diagrams. Before beginning, we emphasize that we discuss Feynman diagrams for 
correlation functions. In particular, following the general point of view common in curved 
spacetime quantum field theory (see e.g. [13]), we view the theory as being defined by 
its gauge-invariant correlators, with the possible existence of a de Sitter S-matrix being a 
secondary issue to be investigated at a later stage. 

Feynman diagrams in Lorentz signature involve integrating products of propagators over 
the relevant spacetime, here dSp. Despite the above exponential decay of de Sitter propa- 
gators, this leads to IR divergences due to the exponential growth of the de Sitter volume 
element ~ (cosht) d . For complimentary series fields with a near zero, even the product of 
3 or more propagators decays only very slowly so that the most familiar tree-level diagrams 
(shown in figure 2) diverge. 

4 



a :- 



d 

"2 + 



- M 



2 „2 



Im 



3" 



> 



-d/2 

i 
i 

V 



Rc 
— 



FIG. 1. On-shell values of a in the complex plane for massive scalar fields. The solid line denotes 
the path of a for increasing M 2 starting from at a = for M 2 = 0. The dotted line shows the 
path of —(a + d) for increasing M 2 starting from — (a + d) = —d for M 2 = 0. Relatively light fields 
with < M 2 t 2 < d 2 /4 correspond to values of a and —(a + d) on the negative real axis and belong 
to the complementary series. Heavier fields with M 2 £ 2 > d 2 /4 correspond to complex values of a 
and —(a + d) on the line Fp and belong to the principal series. 



FIG. 2. Even the above tree-level diagrams diverge for complimentary series fields with a close 
enough to zero. 

Furthermore, even in the principal series, IR divergences arise in generic loop diagrams. 
Consider for simplicity a correction to the propagator; i.e., a diagram with two external 
lines. Let us work in position space and fix the spacetime points x±, xi associated with each 
external line. Then all vertices must be integrated over dSo- Consider in particular the 
integral over the far future region with the relative positions of the vertices held fixed. Then 
the measure contributes a factor of the de Sitter volume e dt but, since the relative positions of 
the vertices are held fixed, the integrand is suppressed only by the propagators corresponding 
to external lines. From the behaviors quoted above, we see that each contributes a factor of 
at best e~ dt l 2 . As a result, the integral diverges at least as J dt ~ t. In particular, so long 

5 



as both external lines describe fields of the same mass, the integrand contains terms that do 
not oscillate at large t. 

It is interesting to note that the above argument applies even to tree-level corrections to 
propagators; i.e., to any quadratic terms in the Lagrangian (mass terms or kinetic energy 
terms) which we choose to treat via perturbation theory. In this context, the above diver- 
gences are related to what was termed a failure of the composition principle by Polyakov 
in [4] - see also [14]. Of course, despite the divergence of the naive Feynman diagrams, 
corrections of this form can always be dealt with by simply diagonalizing the quadratic 
part of the Lagrangian and writing down the resulting free propagators. At least in this 
particular case it is clear that there is no problem with the theory itself, but merely with 
the method of calculation 1 . There is some rough similarity here to the familiar problem of 
secular divergences in classical mechanics, where an infra-red effect appears to be large due 
to not properly accounting for finite shifts in frequency 2 . The relevant question is whether 
another method of calculation can remove all IR divergences and define a useful de Sitter- 
invariant vacuum for the interacting quantum fields. A natural candidate based on analytic 
continuation from Euclidean signature is discussed in section III below. 



III. ANALYTIC CONTINUATIONS IN DE SITTER FIELD THEORY 

It is well known that the Euclidean sphere S D is related to dSp by analytic continuation. 
In particular, the standard metric 

dQ 2 D =f [dr 2 + (sin r) 2 dn 2 d ] , (4) 

on S D (1) can be obtained via the Wick rotation r — >■ t given by 




Because S is compact, no IR divergences can arise through integrals over S . As a result, 
so long as the linearized field theory admits an SO(D + 1) propagator (i.e., so long as 

1 As emphasized in [4] and further explored in [15], the failure of naive Lorentz-signature perturbation theory 
is associated with the fact that there is no adiabatic theorem in de Sitter space. Even a slow change of 
coupling constants in the distant past typically has finite effects at finite times. This phenomenon is in 
turn due to what is effectively a diverging blueshift due to the rapid contraction of de Sitter space in the 
distant past or, what is equivalent, to the spacelikc nature of the past de Sitter boundary I~ (so that 

geodesies enter the future light cone of a given point on I~ only at finite times). 

2 As noted in [16], secular divergences are associateQ even more closely with issues that arise for the special 

case M 2 = which we do not consider here. 



there are no massless scalar fields) , the familiar Euclidean-signature Feynman diagrams will 
converge to define an interacting SO(D + l)-invariant state on the sphere. By this we mean 
that, at each order in perturbation theory, the correlators computed on the sphere will be 
invariant under SO(D + 1) and will satisfy the (Euclidean) Schwinger-Dyson equations. See 
e.g. [17] for other work involving interacting quantum field theory on S D . 

It is therefore natural to attempt to define SO(D, l)-invariant states in the Lorentz- 
signature quantum field theory by analytic continuation from Euclidean signature. Recall 
that a state in quantum field theory can be defined by its correlation functions, so it is the 
correlators upon which the analytic continuation must in fact be performed. Any set of cor- 
relators which satisfies the Schwinger-Dyson equations and appropriate positivity conditions 
may be considered a valid state of the field theory. But by the usual arguments the Lorentz- 
signature Schwinger-Dyson equations are just the analytic continuation of those in Euclidean 
signature, so the Schwinger-Dyson condition is automatically fulfilled by correlators contin- 
ued from Euclidean signature. Similarly, the resulting Lorentz-signature correlators will be 
invariant under SO(D, 1). Furthermore, in the limit of small couplings, positivity conditions 
satisfied for free fields cannot be violated by adding perturbative corrections. In addition, 
the analogue [18] for de Sitter of the Osterwalder-Schrader reconstruction theorem (see e.g. 
[19]) states that this positivity is guaranteed by reflection-positivity of the Euclidean corre- 
lators, which holds at least formally when the potential is bounded below (and which holds 
rigorously for polynomial potentials bounded below if D = 2) [20]. 

In order to gain more intuition for this procedure, it is useful to describe an alternate 
(though computationally more difficult) construction of our state. Because our Lorentz- 
signature correlators satisfy the Lorentz-signature Schwinger-Dyson equations, they may be 
thought of as the result of time-evolving initial data from t — 0. But at t = no analytic 
continuation is required; the Lorentz-signature correlators are precisely the same as that 
Euclidean correlators up to factors of i associated with explicit time derivatives. So our 
Wick-rotated state is identical to what one might call the de Sitter Hartle-Hawking vacuum 
[21] defined by using the Euclidean path integral to compute the state on the S* 3 at t = 
and then evolving away from t — using the equations of motion. 

The existence of perturbative de Sitter-invariant vacuua for interacting (massive) scalar 
field theories is therefore clear at an abstract level. In section IV below, we demonstrate that 
there are no hidden subtleties by computing tree and 1-loop corrections to propagators in 
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precisely this way. Somewhat less trivially, we also explore the large t behavior of the results 
in order to probe the stability and other IR properties of the resulting de Sitter- invariant 
states. Most of our effort will be associated with writing the results in a form appropriate for 
controlling the analytic continuations at large timelike separations. The final output will be 
an integral over functions on dSo which allows us to read off the large t asymptotics of the 
Lorentz-signature correlators. Since the integral converges absolutely, it provides a basis for 
practical numerical calculations even when it cannot be evaluated exactly. In the remainder 
of this section we review two tools that will prove useful in performing the desired analytic 
continuations. 

A. A tool in position space: Embedding distance 

Our focus in this work is on loop-corrected 2-point functions. Since the vacuum on S D 
is invariant under the action of the isometry group SO(D + 1), Euclidean 2-point functions 
{(f>i(xi)<f>j(x2)) may be written as functions of the geodesic distance between xi and x 2 . It 
turns out to be even more convenient to parameterize this separation by the embedding 
distance, i.e. the length of the chord between x 1 and x 2 in an ambient space IR D+1 . This 
embedding distance may be written in terms of coordinates on the sphere as 

Z 12 := Z(xi,x 2 ) = cos T\ cosr 2 + sinri sinr 2 (xi • x 2 ), (6) 

where x\ and x 2 are unit vectors on the sub-sphere S D ~ X . The distance Z is confined to the 
range [—1,1] with 1(— 1) the podal(anti-podal) point. 

Under the analytic continuation (5) the spherical embedding distance (6) becomes the 
SO(D, l)-invariant de Sitter embedding distance 

Zi 2 = — sinhti sinht 2 + coshti cosht 2 (5'i • x 2 ), (7) 

where the embedding space is in this case M 15 ' 1 . On dSp, the values of Z 12 range over all 
of R; in particular, the embedding distance satisfies i) Z 12 e [—1, 1) for spacelike separa- 
tions, ii) Z\ 2 = 1 at coincident points, and iii) |Zi 2 | > 1 for timelike separations. As a 
result, a Euclidean correlation function {(f>i{x\)(f)j{x 2 )) = (0j0 J (Zi 2 )) may be continued to 
the Lorentzian correlator {(pi(xi)(pj(x 2 )) L = (0 i 0j(Zi 2 )) L (or its time-ordered counterpart) 
simply by continuing Z 12 from [—1,1] to IR. 
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Of course, one must deal appropriately with branch cuts and singularities for the result 
to have the desired physical properties. As may be inferred from the flat-space limit, the 
correct definition is 

{T<f>( Xl )(f>(x2)) L := ((j>i<f>i(Z 12 )) , (8) 

where 

Z 12 := Z 12 + ie. (9) 
Similarly, one may define the Lorentz-signature Wightman 2-point function by 

(0(x 1 )0(x 2 )) L :=(^(Z 12 )), (10) 

where 

~ f +ie if X\ is in the future of x 2 

Z 12 = Z 12 +\ . (11) 

I — ie if Xi is in the past of x 2 

B. A tool in momentum space: Watson- Sommerfeld transformations 

As in flat space, de Sitter calculations are typically easiest to perform in what is effectively 
a momentum space representation. Now, it is well known that there are various subtleties 
regarding the definition of de Sitter momentum space in Lorentz signature. For example, 
the spectrum of the wave operator on L 2 (o?S'd) has both continuous and discrete parts 
[12]. However, one may avoid all such issues by simply calculating Feynman diagrams in 
Euclidean signature and using the basis of L 2 (S D ) given by the standard spherical harmonics 
Yg satisfying [22] 

?V 2 x Y £ {x) = -L(L + d)Y £ {x), (12) 

where V 2 is the standard scalar Laplacian on S D with metric (4). Here L — (L D , L D _\, . . . , L ± ) 
is the set of D angular momenta; the Lj are integers satisfying L D > L D ^i > ■ ■ ■ > L 2 > \Li\. 
We will refer to Lp as the total angular momentum. The harmonics satisfy the orthonor- 
mality and completeness relations 

Y,Y t {x)Y* z {y) = £ D ~5(x,y), [ Y z {x)Yt{x) = t D 5 za . (13) 
1 Jx 
In addition, the harmonics satisfy the following very useful relation [23]: 

E = r( ^" gf(^). («) 

3 
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where here L = (L,j), C"{x) is a Gegenbauer polynomial (see appendix A 2), and Z, 



xy — 



Z(x, y) denotes (6) with arguments x, y. 

The usual operations readily express Feynman diagrams on S D as sums over spherical 
harmonics or, equivalently, over Gegenbauer polynomials using (14). One might try to 
obtain useful expressions for de Sitter correlators by analytically continuing such sums over 
polynomials using (9). However, the C£(x) are polynomials for integer L, so each term in 
such a sum diverges at large Z xy . But this is precisely the region we want to study, since we 
wish to determine the behavior of correlators at large timelike separations. 

It is therefore useful to rewrite sums of Gegenbauer polynomials C£ as integrals of more 
general Gegenbauer functions C"(x) over an appropriate contour C in the complex plane 
using a procedure that one might think of as analytic continuation in momentum space. 
Specifically, we use a Watson-Sommerfeld transformation (see e.g. [24]): To express a general 
sum S = J2l s (L) as a contour integral in the complex L plane, one first chooses any function 
s such that i) s agrees s at the values of L appearing in the original sum and ii) s is analytic 
in some open neighborhood of the complex L-plane around each such point. One then 
multiplies s(L) by a meromorphic kernel k(L) having unit-residue poles at the values of L 
appearing in the original sum. One then need only choose an appropriate contour C along 
which to integrate: 



Finally, one may attempt to deform the original contour Co to another contour C over which 
one has more control. 

In our applications, the summand s(L) contains a factor of C^ 2 (x). We therefore take 
s(L) to contain a similar factor in which C^ 2 (x) is a Gegenbauer function (see appendix 
A 2 for conventions) for general complex L. Recall that we wish to evaluate our Feynman 
diagrams at large \Z\. It is therefore useful to know that for general complex L, C^ 2 (Z) is 
a sum of two terms that behave for large real Z like Z L and Z~^ L+d ^ (see appendix A 2). As 
a result, we achieve the most control if we can deform the contour to the line Tp associated 
with principal series values of a (i.e., on which Re L = —d/2) where both terms decay at 
large \Z\ like \Z\~ d / 2 . Our basic goal 3 is to express all diagrams in terms of integrals over 
Tp, and to carefully study the extra terms that arise as one deforms the contour from C 
3 As we will discuss in section IV A, one can obtain even more information about the large \Z\ behavior by 
applying additional tricks, but such embellishments are not needed for the most central results. 




(15) 
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to Tp. If the integrand decays sufficiently rapidly at large \L\, then there is no contribution 
from infinity. The Lorentz-signature propagator will then decay at large values of \Z\ if all 
singularities encountered are sufficiently close to Tp. 

It is useful to quickly illustrate this technique by computing the free propagator. Recall 
that the free propagator A a xy on the sphere is the unique solution to the inhomogeneous 
Klein- Gordon equation 

-(V 2 - M 2 )A° xy = -(V 2 - M 2 )A xy = 5 xy . (16) 

From the above-mentioned properties of spherical harmonics we immediately see that A xy 
may be written 4 

L L 

where in the second equality we've defined 

A Lct := L(L + d) + M 2 £ 2 = (L - a)(L + a + d). (18) 

The expression (17) provides a spectral representation of A° xy on the space (x, y) G S D x S D . 
Other representations may be found by summing over the angular momenta. First, by using 
(14) to sum over all but the total angular momentum one obtains 

^) = ^ D ^i ( ^Oi'\Z). (19) 

This is a spectral representation on the interval Z e [—1,1]- The fact that A a (Z) depends 
only on the invariant distance Z is manifest. Note that in the form (19) one may readily 
extend the definition of A a {Z) to arbitrary real dimensions d, a procedure that will prove 
useful below and for dimensional regularization of UV divergences. 

One now wishes to perform the final sum in (24). To do so, we take s(L) to be 

s{L) = ^±^ e -- L Cf 2 (-Z); (20) 

this is just the summand in (19) rewritten slightly by using the Gegenbauer reflection formula 
(A9). We let 

fc ( i » = £^i) = - £ " tr '- i ' 1 + L '- < 21 > 



4 Dimensional analysis shows the length dimensions [. . . ] of the following quantities: [M 2 ] = —2, [<j)\ = 
^1T^ [-9™] = ni ~ D ^ — D, where g n denotes the coupling constant of n-ficld interactions. It follows that 
[A a ] = 2-D. 
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which inserts poles of unit residue at all L e Z and write 



A CT (Z) 



/2 -D r (I) ("I) 
4^/2+1 r 

rfL (2L + d) 



x 



r [-L, L + d] 2 F 1 



-L, L + d; 



d+1 l + Z 



(22) 



Here 2 Fi(a,b;c; z) is the hypergeometric function and we use a condensed notation for 
gamma functions presented in appendix A 1. Since the hypergeometric function is singular 
at Z = 1, the above procedure should be performed with with Z < 1; we will later continue 
to \Z\ > 1. 

The contour C\ is depicted in figure III B. The integrand has poles at L = 0, 1, 2, ... , 
L = —d, —(d+1), —(d+2), . . . , and at L = a, — (a+d). In fact, the integrand is antisymmetric 
under the reflection L — > —(L + d). Letting \L\ — > oo in the neighborhood of the real axis 
the integrand decays as a negative power of \L\ - this is basically a result of the fact that the 
original series (19) converges - but, since the factors e ±t7rL in (20) and (21) cancel in (22), 
the integrand decays exponentially away from the real axis, i.e. at large |ImL| ^> 1 the 
integrand decays like e _7r l Imi L Because of this, we may deform the contour of integration 
at infinity. Consider deforming the contour C\ to the contour C 2 defined by a straight line 
at any angle to the real axis passing through the reflection point L = —d/2. In the process 
we deform the contour through exactly one of the two poles at L = a or L = —(a + d), 
picking up a residue (see Fig. Ill B). The remaining line integral along C 2 vanishes due to 
the antisymmetry of the integrand under L — > — (L + d). Thus we find that 



A a (Z) = £ 



2-D 



r(|) 



i 



xRes 



4^/2+1 r (d) 
(2L + d) 



A 



La 



r [-L, L + d] 2 F 1 



-L, L + d; 



d+1 l + Z 



(23) 



L=a or — (a+d) 



The residues at L = a and L = — (a + d) are equal (again because the integrand is antisym- 
metric under L — > — (L + d)), so one readily obtains 



A a (Z) = i 



2-D 



1 



-I 



2-D 



4^/2+1 

r(i) 



|, -a, a + d 



-a , a + d; 



d+1 1 



-c d J\-z), 



(24) 



(25) 



47T d / 2 sin(7T(j) 

where the definition of the Gegenbauer function (A3) was used in the final step. In this 
form, the propagator is readily continued to all Z6l. In particular, the large \Z\ behavior 
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FIG. 3. An example of the contour prescription for computing A CT (Z). The contour Ci is an 
arbitrary straight line through the reflection point L = —d/2. Sample cr-poles are drawn for the 
principal series (boxes) and complementary series (circles). 

follows from (A5) which shows that A CT is a sum of two terms, respectively proportional to 

Z a and Z-^ +d \ 



IV. PERTURB ATIVE CORRECTIONS 

We are now ready to compute perturbative corrections to propagators. At tree level, 
such corrections can arise only through interactions of the form £ int = g(f>i(x)(f>2(x), where 
the fields need not be distinct and g has length units [g] = —2. This interaction can be 
thought of as resulting from a non-diagonal mass matrix. At the level of the Lagrangian one 
may of course use a field redefinition to rewrite this theory in terms of free massive scalar 
fields, after which one may compute the corrected propagator exactly in either Euclidean 
or Lorentzian signature. It is nevertheless useful to understand the perturbative treatment 
of such terms, not least because they arise as counterterms needed for renormalization. We 
therefore briefly discuss such corrections before turning to 1-loop corrections in section IV A. 

13 



01 



0"2 



0"1 



<?2 



(a) 



(b) 



FIG. 4. Tree-level contributions to (a) {4>i{xi)4>2{x2)) and (b) {4>i(xi)(f)i(x2))- 
On the sphere one may readily compute the leading correction (figure IVa): 

(&(*i)fc(s 2 )> = -9 f A£A£ + 0(<?) 3 

= M ? - Mf ( A<T1 ( Zl2 ) " ACT2 (^)) + 0{ 9 f. (26) 

The final equality follows from the equation of motion and the fact that no surface terms 
arise upon integrating by parts (the latter statement is not true in de Sitter). The degenerate 
case M\ = Mf (e.g., where <pi and 2 represent the same field) can be found by taking the 
limit of (26); the right hand side becomes —d M 2A a (Z 12 )\ M 2 =M 2. Explicit expressions for 
d a C"(Z) for half-integer and integer a may be found in [25] if desired. The result of course 
agrees with what one finds by diagonalizing the quadratic term in the action, computing the 
exact propagator, and then expanding the result perturbatively in g. 

However, it is useful to note that (26) can also be obtained following the computational 
scheme outlined in section III B. One simply uses (17) and the orthogonality of spherical 
harmonics to convert the integral over S D to a sum over L, and then uses (14) to write 
the the result as a sum over Gegenbauer polynomials. Finally, one may use a Watson- 
Sommerfeld transformation very similar to the one described in section III B to perform the 
final sum over L and obtain (26). In the degenerate case the derivative d a A a (Zi 2 ) arises 
from evaluating a double pole. 

No matter how (26) is obtained, the result is straightforward to continue to Lorentzian 
de Sitter using (9). Again, the result agrees with the correction found by diagonalizing 
the quadratic term in the action, computing the exact propagator, and then expanding the 
result perturbatively in g. This makes it clear that any term of the form (26) is precisely a 
correction to the scalar field mass matrix, and that the corrected propagators will fall-off at 
large Z in the manner one would expect. 
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Q 



(b) 



(a) 



FIG. 5. 0(g) contributions to (0i(xi)0i(x2))- Diagram (a) is the 1-loop contribution while (b) is 
a possible counterterm. 

A. 1-Loop contributions 

We now analyze the possible effects of interactions at 1-loop order. At this level, only 
three- and four-particle interactions can contribute to 2-point functions. For 4-particle 
interactions of the form |(0i(a;)) 2 (02(a;)) 2 (with <f>i and 02 perhaps representing the same 
field), the relevant diagrams are those of figure IV A. Both of these diagrams are of the form 
(26) discussed above. This is manifest for diagram (b), while it becomes clear for diagram 
(a) by writing 



In other words, diagram (a) of figure IV A is just a constant (given by the propagator 
at coincident points) times diagram (a) for figure III B. After renormalizing the constant 
(by using dimensional regularization and perhaps adding the counterterm associated with 
diagram (b) of figure IV A for the case mi = m 2 ), the result just a (real) correction to the 
mass of 0i. Thus the corrected propagators again fall off at large \Z\ as one would expect. 
One can in fact set the mass corrections to zero by an appropriate choice of renormalization 
scheme. Interested readers may find detailed results for the minimal subtraction (MS) 
scheme listed in appendix C for dimensions D = 3,4. 

We therefore turn directly to the more interesting case of 3-particle interactions, which 
provide a more computationally difficult example of quantum corrections to the 2-point 
function. We consider a theory with three massive scalar fields and interactions given by 




(27) 



int 



g<i>i(x) <j) 2 (x)(f) 3 (x) 




(28) 
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01 0-1 

X 

(c) 

FIG. 6. 0(<? 2 ) corrections to (<^i(xi)<^>i(x2)). Diagram (a) is the 1-loop contribution, (b) the 
counterterm due to field renormalization, and (c) the counterterm due to mass renormalization. 
The slash in diagram (b) denotes the action of V 2 . 

As usual, the results for self-interacting fields can be obtained at the end by taking degenerate 
limits where one or more masses coincide. 

The first term in (28) provides the 3-particle interaction while the remaining terms are 
counterterms which arise from the renormalization of the fields and bare masses. As for 
4-particle interactions, we can ignore renormalization of the coupling g as it plays no part 
in the renormalization of the 2-point function at this level. The coefficients in (28) have 
length units [g\ = 2=2, [(Sfa)] = 0, \{5Mf)} = -2. The interaction in (28) is relevant in 
spacetime dimension D < 6 and marginal in D = 6; we will therefore study this theory in 
D = 3,4,5,6. 

The lowest-order corrections to the scalar 2-point function (0i(xi)0i(x2)) occur at 0(g 2 ); 
they are shown diagrammatically in figure IV A. Using the notation of this figure we write 

(M^)M^)) {2) = (a) + (b) + (c). (29) 

Let us first compute diagram (a) on the sphere: 

(a)=g 2 [ [ (30) 

To proceed, we expand the product of two A distributions on S D in the basis given by 
spherical harmonics: 

A-A-(Z 12 ) = f-™J2p^ 2 (L)Y £ ( Xl )Yl(x 2 ) (31) 

L 

= i4 ~ 2D ^Ti 5> + «) Pirin (L)CZ(Z 12 ). (32) 

L=0 

The spectral function p Ul(T2 (L) defined by (31) will be discussed shortly and is computed in 
appendix B. As in this appendix, it is convenient to keep track of the spacetime dimension 
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through the quantity a := d/2 — (D — l)/2. Inserting (31) into (30) and using (17), (13), 
and (14) we find 

r„\ — y8-2D„2 T (q) \^ (L + a)p a2(T3 (L) 



^ (A LCT1 ) 

The counterterms (6) and (c) are straightforward to compute 

J 7/1 



r» ^ (L + a)(^)L(L + 2a) 
^TlZ^ m ^ ^lI^J, (34j 



2tt 



(c) = - (^M, 2 ) I A- A- = 



L=0 

oo 



(Al ct1 ) 2 



D T(a) f^ (L + a) (<5M 2 ) 
Combining our results we have the following expression for (0i(xi)0i(x2))^ 2 ' 

(0i(^)0i(^ 2 )) (2) = ^- D ^ E + 

where we've defined 

9 2 £ 6 ~ D P^ 3 (L) - i 2 mi) - L(L + 2a){8fa) 



(35) 



f(L) := 



(36) 



(37) 



Let us now discuss the function p aiU2 (L) defined by the expansion (31). From the orthog- 
onality of Gegenbauer polynomials we may compute 

27T Q+1 



p ai a 2 (L) := t 



2D-4 



l - [ +1 dZ(l-Z 2 ) a - 1 / 2 Ct(Z)A^(Z)A^(Z), (38) 

l L J-l 



T (a) (L + a) Aj 

where is the Gegenbauer normalization (A8). The function p aia2 (L) is clearly invariant 
under the actions 

<7i ->■ -((7i + 2a), cr 2 ->■ -(cr 2 + 2a), <7i < — >■ cr 2 . (39) 

Near Z = 1 the distribution A CT (Z) behaves like ~ (1 — Z) 1 / 2 ^, so we see that the integral 
in (38) converges for < a < 3/2. We compute p aia2 (L) for this range of a in appendix B; 
the result may be written 

, 1 COs(7T<7i' 



x r 



x 7 V 6 



I6n a sin 71(0"! + a) 

2 - 2a, -(7i, L + 1, 2 + L — <7i — a, L ~^ 2 , L '^+^+ 2 " 
1 — <7i — a, L + a + 1, L + 1 — <7i, 



2 ' 2 

4+L— cri — g-2— 4a 4+L— cr\+rJ2 — 2a 



1 + L — o"! — a; 1 — a, 1 — <7i — 2a, 1 + L, 



2 ' 2 

L — o"i — c 2 L — cri + <7 2 + 2a 



+(<Ti ->■ -(<Ti + 2a)) ^ + (<7i < — ► a 2 ) 



(40) 
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Here 7Vq{cl; b, c, d, e, f) is a so-called very well-poised jFq hypergeometric function (see ap- 
pendix A 3). The series defining the -jV^a; b, c, d, e, /) in (40) is absolutely convergent for all 
complex L, a, a±, and <r 2 . We may define p aia2 (L) for complex a via the analytic continua- 
tion of (40) beyond the interval < a < 3/2. As is discussed in appendix B3, this extended 
p aiU2 (L) has poles at a = 3/2,5/2, .... We also show in this appendix that p aiU2 (L) has 
poles in the complex L-plane at 

L = o\ + 0"2 — 2n, —o\ + o" 2 — 2a — 2n, o\ — a 2 — 2a — 2n, —U\ — a 2 — 4a — 2n. (41) 

for n G N . We will address the meaning of these poles momentarily. An important property 
of p aia2 (L) defined in (40) is that it obeys 



p ai a 2 (L) = pcTa^(L) = p ai a 2 (L) 



(42) 



for "on-shell" masses U\ and <j 2 . The first equality follows from the fact p aia2 (L) can be 
written as an absolutely convergent series of terms which may be expressed in terms of 



Gamma functions, and the Gamma function itself obeys T(x) = T(x). The second equality 
follows for on-shell values of <7i and er 2 . "On-shell" values of a are either (i.) —a < a < 0, 
in which case a G R, or (ii.) a = —a + iv, v G R, for which a = —a — iv — —{a + 2a). 

We can now discuss the renormalization coefficients in (37). We use these coefficients to 
cancel any superficial divergences in p (T2(T3 (L) and render f(L) finite. For the dimensions of 
interest, such superficial divergences occur when a = 3/2 and a = 5/2 (D — 4 and D = 6). 
In the neighborhood a = (3 — e)/2, p a2a3 (L) diverges as 

1 



(L) 



87T6 



+ 0(6°). 



(43) 



a=(3-e)/2 

Following the MS scheme, this divergence is cancelled by setting 

-.2 



a={3-e)/2 

For a = (5 — e)/2 we have 



9 

87T6 



+ 0(g*), (6<f>i) 



(44) 



a=(3-e)/2 



0"2 <T3 



(L) 



1 



647r 3 e 



L(L + 5) 
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a=(5-e)/2 

This divergence is cancelled by setting 
= _S 2 (Mf + Mf 

a=(5-e)/2 



+ 0(e°). 



(45) 



647r 3 e 



+ 0(^ 4 ), (50.) 



a=(5-e)/2 



1927T 3 e 



+ 0(^ 4 ). (46) 
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The expressions for (<5M|) and (<5Mf ) are given by the obvious permutation of the masses. 
For D = 3, 5 we set (5Mf) = (Sfa) = 0(g 4 ). 

Let us now return to our expression (36) for {4>i(xi)(f)i(x2)) (y2 \ Our task is to rewrite 
this in a form well-suited to analytic continuation to Lorentz signature. We proceed in 
the same way we dealt with the free 2-point function in (19), using a Sommerfeld- Watson 
transformation defined by the same kernel (21) and integrating along a contour C enclosing 
the poles at L — 0, 1, 2, 

(M^)M^)) {2) = -2<f ^-.f(L)(L + a)A\Z 12 ), (47) 

where we have used (24) to replace Gegenbauer functions by A distributions for general real 
a = d/2. 

The integrand decays exponentially away from the imaginary axis like e _7r l ImL l; we can 
therefore deform the integration contour away from C . We would like to deform the integra- 
tion contour to the contour V along the straight line Tp defined by Re(L) = —a (see figure 
IV A). By convention, we take Y to pass on the left side of any poles that lie precisely on 
T P . 

As we deform the contour, we will pick up residues from any poles we encounter. The 
integrand in (47) has many poles in the L-plane. The distribution A L (Z) has simple poles 

at 

L = n, L = -(n + 2a), for n G N . (48) 

The function f(L) has the simple poles in p - 2cr3 (L) listed in (41); in addition, the (Ai CT1 ) 2 in 
the denominator of f(L) has double-poles at 5 

L = <7i, L = -(<ti + 2a). (49) 

Despite all these poles, only a very few poles are encountered as we move the integration 
contour from C to Y. When <f>i(x) is in the complementary series then —a < U\ < and the 
pole at L — o\ is on the right-hand side of T. When 4>\(x) is in the principal series both the 
poles at both L = o\ and L = — (<7i + 2a) lie on the line Tp. Additionally, if both faix) and 
4>s(x) are in the complementary series it may be that —a < o<i + 03 < and even possibly 
5 Furthermore, there are the special cases where u\ = — (a\ + 2a) = —a; in this case, due to the (L + a) in 
the numerator of the integrand, there is a 3rd-order pole at this point. There is also the possibility that 
01 = a 2 + 03 1 in which case a double-pole exists at this point. Both of these special cases can be found 
as limiting cases of the more general case so we will not treat them explicitly. 
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—a < a 2 + a 3 — 2 < 0; in these cases the poles at L = a 2 + cr 3 and L = (T 2 + o" 3 = 2 lie to 
the right-hand side of T. We conclude that 



<0i(-i)0i(* 2 )> (2) = 2Res{/(L)^^ 

+2 J^J(L)(L + a)A L (Z 12 ). (50) 



Here an asterisk notes that the residue should only be considered if the pole location has 
Re(L) > —a. See figure IV A for examples. 

Let us first consider the case when <pi(x) is a complementary series field. In this case 
—a < <Ji < 0. It is straightforward to evaluate the residue at L = ax, which may be written 
in the useful form 



R, :=2Res{f(L)(L + a)A L (Z 12 )} L=ai 
d 



d{M 2 ) 



[g 2 Re {i 4 ~ D p^MM 2 ))} - (6M*) + M 2 (^)] ^ {M ) {Z l2 ) 



M 2 =Ml 



(.51) 




FIG. 7. Two examples of the complex L-plane. Poles at o\ and — (ai+2a) are marked by CPs, while 
the pole at (a 2 + 03) is marked by an o. In example (a) 4>\{x) and (fi 2 (x) are in the complementary 
series while (f>3(x) is in the principal series. Here r = Tp. In example (b) (ft 2 (x) and 03 (x) are in 
the complementary series while (j>i(x) is in the principal series. Here T differs from Tp since we 
take T to pass to the left of poles with Re L = —a. 



20 



Next we examine the integral over T. Inserting L = —a + iv we have 

I:=2 j r ^ + «)/W AL (^12) 
• roc 

= dvv [f(-a + iv)A' a+iv {Z l2 ) - f(-a - iv) A~^{Z X2 )\ 

71 Jo 



ig 



2p6-D 



oo 



dv 



2g 



^ Jo 

2 fi2-D :■ x 



(A- 



) 2 



dv 



v lm{p a2as (-a + iv)} A _ a+il/ 



M 2 ■ ) 2 

lvl -a+iv) 



(Zu). 



(52) 



Jo {Ml 

The first equality merely uses the symmetry of the contour under complex conjugation to 
write the expression as the integral of a quantity that is manifestly real (for real Zi 2 ). The 
second equality then follows by using the relations A_ Q _j I/(T1 = A_ Q , +iiyiCri and A~ a ~ zu (Z) = 
A- a+iu {Z), the definition of f(L), and the property (42). 

Let us now consider the case when (f>i(x) is in the principal series. In this case o\ = 
—a + ivi, v\ G R, so both the poles at L = o\ and L = — (<Ji + 2a) lie along the line Tp. The 
final result is almost identical to the result for the case of complementary 0i(x). Recall that 
the contour T is indented as shown in figure IV A (b). Thus we have two poles whose residues 
combine to give the twice the expression on the right-hand side of (51). Furthermore, the 
integral over T is of the same form as / (52), but with the contour deformed slightly to the 
left. It is convenient to remove the indentations in the contour by instead writing the result 
as the principal part of / added to (51), where the deformation of T back to Tp precisely 
compensates for the extra factor of 2 noted above. 

Finally, we must consider the case where both <p 2 (x) and 3 (x) are complementary series 
fields with sufficiently light masses such that —a < cr 2 + cr 3 < and possibly —a < a 2 + 
o"3 — 2 < 0. In these cases we encounter pole(s) at L = a 2 + 03 (and L = a 2 + 03 — 2) as we 
move the contour. These residues are easily evaluated using (B26): 



R 2 :=2Res{f(L)(L + a)A L (Z 12 )} 
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4tt«+i (A CT2+CT3 , CT1 )2 
R 3 := 2Res{f(L)(L + a)A L (Z 12 )} L=a2+aa _ 2 



L=(72+<T3 

— cr 2 , 0"2 + — 0"3> °3 + Ot 

—o 2 — (T 3 , a 2 + 03 + a 



A CT2+(T3 (Z 12 ), 



(53) 



2g6-D 



a(a 2 + 0-3 + 2a - 2) 



vr Q+1 (A CT2+a3 _ 2i(J1 ) 2 (<7 2 + a 3 + a - 1) 



1 - cr 2 , <7 2 + a - 1, 1 - <7 3 , 03 + a - 1 

2 — a 2 — cr 3 , <r 2 + a 3 + a — 2 



(54) 
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Assembling our results we have the final expression 

(0i(zi)<Mx 2 )> (2) =Ri + P(I) + R 2 + R 3 , 



(55) 



where Ri, I, R 2 , R3 are given respectively in (51), (52), (53), and (54). In (55) the R 2 
term should be included only when —a < a 2 + 03 and likewise the R 3 term should only be 
included when —a < o 2 + 03 — 2. The P in (55) is a reminder to take the principal part in 
integrating through any pole terms on the axis in the integral I. This result is manifestly 
real for on-shell masses as it should be. Earlier we noted in a footnote that that there are 
two degenerate cases in which the computation above requires modification, namely when 
G\ = —a and when G\ = o 2 + 03. One can find the correct result for these cases by taking 

(2) 

the appropriate limits of (55). Finally, the Lorentz-signature correlator (T(f) 1 {xi)(f)i{X2))\ j is 
defined as (55) with Z\ 2 — > Z 12 ; likewise, we define ((f>i(xi)(f)i(x 2 ))^ as (55) with Z 12 — > Z 12 . 

Our final expression (55) is rather complicated. However, it has two very useful features. 
The first is that the remaining integral / converges absolutely for arbitrary \Z 12 \ > 1 so 
long as the contour has been deformed away from any poles in the appropriate manner to 
compute the principal part. As such, the result is amenable to numerical calculations and 
gives a practical tool for extracting detailed physics. The second is that it allows us to 
extract the large \Zi 2 \ ^> 1 behavior and so to study the corrected propagator in the deep 
IR. At large \Z 12 \ 3> 1 the first three terms in (55) have leading behavior, in order, 

\Z 12 r\ogZ 12 , \Z 12 \~ a± ™, \Z l2 \^. (56) 

The term with the slowest decay provides the leading behavior at large time-like separation. 

As a slight aside we mention that our final expression (55) can easily be brought into the 
Lehmann-Kallen form of the 2-point function [26] 

poo 

(MxiWiM)^ = / dM 2 p(M 2 )A M2 (Z 12 ), (57) 

where p(M 2 ) (which should not be confused with p aiU2 (L)) is the spectral density, and the 
integral is over positive real M 2 > 0. The integral / in (52) is already in the form of an 
integral over the principal series masses, so one need only encorporate the poles R\ (and 
possibly R 2 and R3) into the integral over M 2 by using delta functions in the obvious manner. 

It is useful to check our results by taking the flat-space limit of (Tfaixjfafa))'?. For 
convenience, let us suppose that no massless fields arise in this limit. Thus all three fields 
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must be in the principal series. For principal series fields we have only the terms R\ + P(I). 
Now, since A a (Z 12 ) reduces in this limit to the flat-space propagator D m2 (xi — x 2 ), we must 
have 

^- D p^ 3 (L)^p^l Mi (k 2 ), (58) 

where L 2 /i 2 — > k 2 and p^\ M 2(k 2 ) is the analogous spectral function of the product of two 
Minkowski propagators defined by 

d mi ( Xi - x 2 )D Mg ( X1 - X2 ) = J ^P^M?s k2 y k{xi ' X2) - ( 59 ) 

One may explicitly check this result for a = 1 and a = 2 where we have simplified ex- 
pressions for p a2(T3 (L). (We have also explicitly verified that the singular 0(1/ e) terms are 
equivalent for a = 3/2 and a = 5/2.) After changing the integration variable in I to 
m 2 = v 2 /£ 2 , one may take the limit £ — » oo holding fixed g, M 2 , (SMf), and (6<f>i). Noting 
that Im{p^2 M |(— m 2 )} = for m 2 < (M 2 + M 3 ) 2 , one finds that the answer agrees with the 
known flat-space result [26, 27]. 

We conclude this section with a brief discussion of correlator given by the sum of single 
particle-irreducible (1PI) Feynman diagrams. This 1PI correlator may be written 6 

<A<MZ 12 )> 1PI = ± GE(z u ), (60) 

where for the diagrams of figure IV A the dimensionless self-energy is 

n(L) := g 2 i 6 ~ 2D P(T2a ,(L) - £ 2 (5M 2 ) - L(L + 2a) (Sfa) . (61) 

This correlator may be analytically continued to de Sitter in essentially the same way as 
the 0(g 2 ) correlator above. An interesting feature we wish to point out is that when 0i 
belongs to the principal series and Re((x 2 + <r 3 ) < —a, the Lorentz-signature (time-ordered 
or Wightman) 1PI correlator decays at large \Zi 2 \ ^> 1 more rapidly than any free 2-point 
function. 

To see this it is convenient to rearrange the summand of (60) slightly before performing 
the Watson-Sommerfeld transformation. We use the Gegenbauer recursion relation (A6) to 
write 

(*MZn)) in = t { - XMm ! m + 2) }cr\z a) , ,62) 

6 Implicit in the sum over 1PI diagrams is the assumption that \H(L)/Mf \ < 1. 
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(<hMZ*)) m - ~^ dT ^ I ^ {x Lai - U (L) - A L+2 , CT1 -U(L + 2) } 



where we also use the fact that C" 2 (Z) = C\[Z) = for the values of a of interest. Using 
our standard Watson-Sommerfeld kernel we may write the sum (62) as 

,r(o) /• ill \ i l | C2 +1 (-z 12 ) 

sin 7i L ' 
(63) 

where 7 is a contour parallel to and slightly to the left of the imaginary axis. Note that at 
large \Z± 2 \ ^> 1 the Gegenbauer function C c ^ +1 {Zi 2 ) behaves like (Z\ 2 ) L and (Zi 2 )~ L ~ 2a ~ 2 as 
compared to C£(Z 12 ) which behaves like (Z 12 ) L and (Z 12 ) _L ~ 2a . The advantage of changing 
the underlying Gegenbauer function from C"(Z 12 ) to C2 +1 (Z 12 ) is that we can shift the 
contour 7 as far to the left as ReL = — (a + 1) while still increasing the decay of the 
integrand at large \Zi 2 \. 

To determine the behavior of the Lorentz-signature correlator at large \Z\ 2 \ ^> 1 it is 
sufficient to determine the first pole encountered as we shift the integration contour 7 to the 
left. For the case of interest the first poles encountered arise from the first term in brackets 
and are located near the on-shell poles L rs o\ and L —o\ — la. Letting o\ = —a + iv 
with v > we may write the on-shell poles in the free theory as L± := —a ± iv. The 
location of the on-shell poles in the interacting theory is given by solving for the zeros of 
the denominator 

X La - n(L) = L(L + 2a) + M\i 2 - U(L) = 0. (64) 

Having computed I1(L) to 0(g 2 ) we may easily solve for the 0(g 2 ) corrections to L±; the 
result is 

Rell(-a + iv) 



ImII(- a + iv) , . 

j± = —a H ± 1 

2v 



v — 



(65) 



2v 

(recall that n(L) is 0(g 2 )). In writing this expression we've made use of the fact, introduced 



above, that p a2(T3 (L) = p a2Cr3 (L). These poles contribute residues to the 1PI correlator 

proportional to C2^ 1 (—Z 12 ); as a result, the leading behavior the the Lorentz-signature 
correlators are given by 

|— Q+Im(— a+iu)/2v±iui (66) 



'12 



with u> = v — Re Pi(—a + iv)/2v. 

The interesting feature of this result is that when Imil(— a+iv) < the Lorentz-signature 
correlator decays faster than \Zi 2 \~ a , i.e. faster than any free 2-point function. The only 
term in Il(— a + iv) that can become imaginary is the function p a2(T3 (—a + iv). For the cases 
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of a — 1 and a = 2 (spacetime dimension D = 3 and D = 5) one may explicitly examine 
the imaginary part of p a2US (—a + iv) using the formulas (B22) and (B23) and find that 
Im(p - 2CT3 (— a + iu)) < 0. Furthermore, in the flat-space limit (large a±, a 2 , cr 3 up to an overall 
scaling of II to this order) it is known that Imll < [26, 27]. For a = 3/2 and a = 5/2 
(D — 4 and .0 = 6) with small <7jza we have performed only a small numerical sampling of 
Im(p - 20 - 3 (— a + iv)), but in all cases have likewise found that Im(p (72(J3 (— a + iv)) < 0. 

We interpret this result as the appearance of a decay channel that occurs when Re(o"2 + 
<7 3 ) < —a and <f>i is in the principal series. Note that the requirement Re(o" 2 + o" 3 ) < —a does 
not in general fix the relative size of Mf + Mf and Mf . Indeed the appearance of this decay 
channel is quite generic for any principal series field <f>i\ only when the intermediate states 
are very light, i.e. satisfying —a < o-i + 03 < 0, does the behavior of the 1PI correlator differ 
from (66). The ability of particles to decay into daughter particles with lighter or heavier 
masses is a natural phenomena in de Sitter space [8-11]. Due to the lack of a globally time- 
like Killing vector field, there is no positive definite conserved energy which would preclude 
such a process. 

V. DISCUSSION 

We have computed loop corrections to Lorentz-signature propagators for de Sitter- 
invariant vacuua in scalar field theories by analytically continuing results from Euclidean 
signature. Our results apply to all masses for which the free Euclidean vacuum is well- 
defined, including values in both the complimentary series and the principal series of 
SO(D, 1). We have provided explicit results in dimensions D > 3 for which the above in- 
teractions are renormalizeable. Our results generally take the form of absolutely convergent 
integral representations sufficient to extract the leading behavior of the Lorentz-signature 
2-point functions at large separations. The absolute convergence implies that such repre- 
sentations are amenable to numerical calculations, demonstrating that our methods provide 
practical tools for calculating Lorentz-signature correlation functions. We have provided a 
number of checks on our results, including consistency with known flat-space limits. Our 
basic methods appear to apply to higher loops as well. 

Of course, our use of perturbation theory requires small couplings. As described recently 
in [28], perturbative corrections in de Sitter space are controlled by a combination of the 
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coupling and the particle masses which diverges in the limit Mi — > 0. In this limit (tak- 
ing all masses equal), we indeed find that the contributions from the 1-loop diagrams are 
proportional to g 4 /(M 4 i D ) or gf/(M 6 £ D ), where g 3i4 are the 3- and 4-point coupling in the 
Lagrangian. Our results for the 4-particle interaction agree with those of [28] before the 
application of dynamical renormalization group (DRG) techniques, though it was shown in 
[28] that DRG resummation can ameliorate the M — > growth to some extent. It would be 
interesting to combine DRG techniques with our Euclidean approach. 

With this caveat, we find that the corrected propagators fall off at large separations at 
least as fast as one would naively expect. Such results are in qualitative agreement with those 
obtained using stochastic inflation techniques [7] , which are expected to be valid in the limit 
Mi <C 1, where i is the de Sitter length scale. Interestingly, for one-loop corrections from 
3-particle interactions we found that, in some cases, the corrected (lPI-summed) propagator 
decays faster than any free propagator with M 2 > 0. This indicates that the vacuum state 
constructed by analytic continuation of all Euclidean correlators is well-behaved in the IR. 
In particular, similar fall-off of higher n-point connected correlators would indicate that this 
state is stable in the following sense: Consider a state \m) constructed by acting on the 
vacuum with (integrals of) m field operators <pi(x) at or near some initial time t = 0. Then 
the n-point functions of \m) are just (integrals of) 2m + n-point functions in our vacuum. 
Let us now consider a limit in which the arguments of such an n-point function retain fixed 
relative separations, but in which each argument is taken to some large time in the future; 
i.e., so that the the n-points are far from the m operators originally used to construct the 
state \m) (which remain near t = 0). Decay of connected correlators means precisely that 
correlators factorize at large separation. Thus, at large times t any n-point function of 
\m) would approach the product of (m\m) with the corresponding n-point function in the 
vacuum. One might say that, when viewed as a functional on local products of quantum 
fields, the large time limit of any above state coincides with our vacuum state. It is natural 
to refer to any such vacuum as being stable. More specifically, when all correlators in a 
given state factorize in the above limit we will say that the state is an attractor state for 
local correlators. 

Strictly speaking, a 3-point function or higher is needed to test this notion of stability, 
while we have computed only propagators here. We will provide a detailed discussion of 
higher n-point functions elsewhere [6], but for now we merely note that our propagator 
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calculations suggest that the perturbative vacuum theory is well-behaved in the IR. 

Supposing that the higher 2m + n-point functions continue to indicate stability of de 
Sitter-invariant vacuua for (massive) scalar field theories, one may be tempted to ask why 
such vacuua should be stable. For familiar vacuua in flat Minkowski space there is a simple 
physical answer: each such vacuum minimizes a positive-definite energy. But de Sitter space 
has no positive-definite conserved energy due to its lack of a globally-defined timelike Killing 
field, so we must search elsewhere for an explanation. The best answer is probably that 
de Sitter space does admit Killing fields that are timelike in a globally hyperbolic domain 
known as the "static patch" associated with that Killing field. Such domains may be treated 
as spacetimes in their own right, with no need to impose extra boundary conditions. For 
positive potentials, the associated Hamiltonian is bounded below in this restricted spacetime. 
As a result, positivity and conservation of this energy forbids instabilities of scalar quantum 
field theories (with positive potentials) in any static patch. Any possible instability of de 
Sitter space must therefore be a more subtle sort, and would not be directly visible to single 
any freely falling observer. 

In any discussion of de Sitter space, it is tempting to ask about quantum gravity effects. 
Although gravitons are massless, they admit a free Euclidean vacuum state [29]. It is 
therefore plausible that our results may generalize to graviton n-point functions (though 
there is a certain tension with the results of [30-33]). Such a result would again preclude 
perturbative instabilities in this context - at least in cases where they are not already present 
at the classical level. However, even in this case there may still be room for more subtle 
quantum gravity effects associated with large regions of de Sitter space (see e.g. [33-37]) 
which are not instabilities per se, and which remain to be investigated in more detail. In 
addition, there are clearly interesting quantum effects involving gravity coupled to scalars 
with very flat potentials. This exception is allowed due to the fact that free massless scalars 
are already marginally unstable (at both the quantum and classical levels). The prime 
example of such an interesting quantum effect is of course eternal inflation, which will occur 
barring the discovery of further novel phenomena, and which may have further implications 
for understanding quantum de Sitter space [35]. 
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Appendix A: Notation and conventions 
1. The Gamma and related function 



We use the following notation in the main text. The Euler Gamma function is denoted 
T (z), and we use the condensed notation 

ai, a 2 , • • • , dj 
h, b 2 , ■ ■ ■ , h 

We also define the Pochhammer symbol for complex a and n G No 

= (a)(a + l)---(a + n-l); (A2) 

thus, (a) n is simply a polynomial of a of order n. The digamma function ip (z) is the 
logarithmic derivative of the Gamma function ip (z) := fh)- 



r( ai )r(a 2 ) •••!>,•) 
r(6!)r(6 2 )...r(6 fc )" 



(Al) 



a + n 
a 



2. Gegenbauer functions and polynomials 



The Gegenbauer function of the first kind may be defined via the hypergeometric function 



Cf(z) := r 



2a + A 
1 + A, 2a 



2^1 



-A , A + 2a ; a + 



1 l-z 



(A3) 



Here a, A, and z are arbitrary complex numbers. Important features of this function, 
including its analytic properties, recursion relations, asymptotic forms, etc., are presented 
in [38]. The function's relation to representations of SO(n) and related groups is nicely 
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described in [12]. Here we present only information used in the text. For \z\ > 1 the 
Gegenbauer function may be usefully rewritten [38] 



C a x {z) 



+ r 



A + 2a, -(A + a) 
a, —A, A + 1 

A + a 
a, A + 1 



A + 2a 1 + A + 2a 



; A + a + 1 ; z 



(2z) x 2 F 1 



A 1 - A 



; —A — a + 1 ; z 



-2 



(A4) 



+ r 



( 2z y(X + 2a) J 1 + (Z- 2 )] 

(2z) x [l + 0(z- 2 )}. 



(A5) 



2 ' 2 

From this we see that at large \z\ ^> 1 the Gegenbauer functions have two asymptotic 
branches, namely, 

A + 2a, -(A + a) 
a, —A, A + 1 

A + a 
a, A + 1 

Gegenbauer functions satisfy many recurrence relations; some that we will make use of are 

(A6) 

d -Cl{z) = 2 n (a) n C^:(z). (A7) 

When A = L G N the hypergeometric series terminates and Gegenbauer functions re- 
duce to the Gegenbauer polynomials. The Gegenbauer polynomials C£(z) form a complete 
orthogonal basis on the interval z E [—1, 1] with respect to the measure (1 — z 2 ) 01 ^ 1 ^ 2 . They 
have normalization 

L + 2a 



(A + a)Ct{z) = a [C« + \z) - C^(z)} , 



1 dz(i - z 2 Y^ 2 ci(z)cUz) = r 

_i {L + a) 



L + 1, a, a 



>LM- 



Gegenbauer polynomials obey the reflection formula 

Cftz) = {-l) L Ct{-z). 
The integral of three Gegenbauer polynomials with common degree a is [12]: 
D(a;L, M, N) := £\z(l - z^-^^C^C^z) 



(A8) 



(A9) 



)l-2a 



7T 



, (A10) 



J + 2a, J - L + a, J - M + a, J - N + a 
J + a + 1, J-L + l, J-M + l, J-N + l 
when J := (L + M + N)/2 G No, and L, M, and iV satisfy the triangle inequalities; otherwise 
D(a,L,M,N) = 0. 
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3. The function 7 14 (a; b, c, d, e, /) 



The function 7 Ve(a; b, c, d, e, f) is an 7 F 6 hypergeometric function with unit argument and 
a special form of the parameters [39] : 



7 V 6 (a; b, c, d, e, /) 

a, 1 + b, 



7-^6 



C, 



d, 



e, 



; 1 



§, 1 + a — 6, 1 + a — c, 1 + a — d, 1 + a — e, 1 + a — /, 



(All) 



The series defining (All) converges when it's parametric "excess" s = 4 + 4a — 2(6 + c + 
d + e + /) has a real part that is greater than zero. The series terminates when one of the 
parameters is a negative integer. When the series terminates because one of b, c, d, e, f is a 
negative integer and the excess takes the value s = 2 the series may be summed and the 
result is known as Dougall's formula: 

t 7- / 7 7 v {l + a) n (l + a-b-c) n (l + a-c-d) n (l + a-b-d) n 

7 V 6 (a; b, c, d, e, -n) = — — — — — — A12 

(1 + a - b) n (l + a - c) n (l + a - d) n (l +a — b — c — d) n 

with e = l + 2a — 6 — c — d + n and n G No- 

There exist a large number of relations between functions of the form -jV^a] b, c, d, e, /) 
different parameters. One such relation of which we will make use is 

7 V 6 (a; 6, c, d, e, /) 

1 + a-e, 1 + a- f, 2 + 2a-b-c-d, 2 + 2a-b-c-d-e-f 
1 + a, 1 + a-e- f, 2 + 2a-b-c-d-e, 2 + 2a-b-c-d-f 
x 7 V 6 (l + 2a-b-c-d; 1 + a-c-d, 1 + a - b - d, 1 + a - b - c, e, /). (A13) 

This equality is valid so long as the series on both sides converge, i.e. that the excess of 
both series is greater than zero. 

It is convenient to define the regularized function 



7 V 6 (a; b,c,d,e, f) 



7 V 6 (a;b,c,d,e, f) 



r[§, 1 + a - b, 1 + a-c, 1 + a-d, 1 + a- e, 1 + a- /] 
(a) n (l + a/2) n (b) n (c) n (d) n (e) n (f)n 



E 

n=0 



r [1 + n, a/2 + n, 1 + a-b + n, ... , 1 + a- f + n] 



(A14) 
(A15) 
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This series defines an entire function in all of it's parameters. Like 7V&{a] b, c, d, e, f) the 
series terminates when one of the parameters is a negative integer. When —a e N , i.e. 
when a = 0, —1, —2, . . . , the series is zero. 



Appendix B: Calculation of A ai A a2 (Z) 

In this appendix we compute the spectral representation (31) of the product of two free 
Euclidean 2-point functions on the sphere S D . As discussed in the main text, this amounts 
to computing (38), where a := d/2. Using the definition of the constant from (A8), one 
may check that p aiCT2 (L) = —p aia2 (—L — 2a) for L e N , 01,02 £ C In this appendix we 
consider only such positive integer L unless otherwise noted. 

From (38) it is clear that p aia , 2 (L) will not in general be finite. Recall that near Z = 1 the 
2-point function diverges as A CT (Z) ~ (1 — Z) 1 l 2 ~ a ) so the integrand (38) diverges near the 
boundary Z — > 1 for a > 3/2. We handle this divergence using dimensional regularization; 
i.e., we consider p aicr2 (L) as a function of the real parameter a, evaluate p ai(T2 (L) for a < 3/2 
for which the integral (38) converges, and then define p Uia2 (L) for a > 3/2 via analytic 
continuation of our final expression. The remainder of this appendix is concerned with 
computing p airT2 (L) and then presenting a number of checks of our work. 

We now turn to evaluating (38) for p aia2 (L). We begin by defining 

2(L + a) _ 2(L + a) _ 1 1 

La ' \ La (L - a)(L + a + 2a) L-a L + a + 2a' { ' 

and inserting (19) twice into (38). We find 

27r a+1 1 T 2 (a) 



T (a) (L + a) A a L (47r«+ 1 ) 2 
00 00 „_|_i 

x 

M=0 N=0 



J2 E A m CTi A^ 2 / dZ{\- Z 2 r^ 2 C a L (Z)C^(Z)CUZ) 

M=0N=0 

r ( ) 1 00 00 
S^(L + a) Ai SI E Am„a^c(«; l, M, N) 
y ' / l M=0N=0 

S aia2 (L). (B2) 



1 



L + l 

a, L + 2a 



To get to the second line we perform the integral using (A10), and in the third line we've 



31 



defined 



M,N 



J + 2a, J - L + a, J - M + a, J - N + a 
J + a + 1, J-L + l, J-M + l, J-N + l 



where, as in (A10), the sum is over all M and N is such that 
L + M + N 



J := 



We can incorporate these restrictions by a change of variables: 

_ -L + M + N T ^ L + M-N Ar 
G : = = J - L, K := — J — N, 



such that 

M = G + K, N = G + L — K, J = G + L. 
In terms of these variables S aia2 (L) becomes 



oo L 



^1(72 (L) — ^2 ^-G+K,<nh-G+L-K,<T2 T 



G=0 K=0 



(B3) 



G N , \L - M\ < N < L + M, \L-N\<M<L + N. (B4) 



(B5) 



(B6) 



K + a, L — K + a, G + a, G + L + 2a 
K + l, L-K + l, G + l, G + L + a + 1 



(B7) 



In the next two sections we sum over first K and then G. 



1. The K-sum 



Let us perform the sum 

L 

H(L; G) := ^ A G+Ki(Tl A G+L _x, CT2 T 



K=0 



K + a, L-K + a 
K + l, L-K + a 



(B8) 



by recasting it as a contour integral in the complex fT-plane. We do so by multiplying the 
summand by 

7rcot(7rX) = - cos(n K)T[-K, K + l]. (B9) 

This function has poles with unit residue at K = 0, 1, . . .,L. Consider then the contour 
integral 

T[K + a, L — K + a] 



I ■■= (-1) j ^— cos(irK)A G+Kj<T1 A G+L _ Kjl 



o. 



L+l 



(BIO) 
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The contour is chosen to be an arc near infinity; because the integrand behaves at large 
\K\ 3> 1 like ~ \K\ 2a ~ 4: the integral vanishes for a < 3/2. By Cauchy's formula it follows 
that the sum of the residues of the poles enclosed in C must likewise sum to zero. The 
integrand has the following simple poles: 

1. K — 0, 1, . . . , L — 1, L, due to (— K) l + \ in the denominator, 

2. K = —a — n, n <E N , due to T (K + a) in the numerator, 

3. K = L + a + n, n & N , due to T (L — K + a) in the numerator, 

4. K — —G + <7i, K — —G — <j\ — 2a, due to Aq+k,^, 

5. K = G + L — (T 2 , K = G + L + a 2 + 2a, due to A g +l-a> 2 • 

We assume for simplicity that <t; ^ —a + Z such that none of the above-mentioned poles 
overlap. There is nothing peculiar about these configurations and we will see that our final 
result is perfectly regular at these values of the a,. 

Let us now turn to evaluating the residues of these poles. 

1. Poles at K — 0, 1, . . . , L — 1, L: By construction the residue of these poles sum to 



H{L-G). 



2. Poles at K = —a — n: These poles sum to the infinite series 



n + a, n + L + 2a 




n + 1, n + L + a + 1 



(Bll) 



n=0 



3. Poles at K = L + a + n: These poles sum to the infinite series 



■oo 



n + a, n + L + 2a 
n + 1, n + L + a + 1 




(B12) 



n=0 



4. Poles at K = —G + <j\, K — —G — a 1 — 2a: These give two terms, 



sin7r(<7i + a) 



7T COS7TO"! 



A-2G+L-<ti,o-2 r 



G — <j\ , G + L — o\ + a 

G + L + \-ox,G + \-ox-a 



+ ((Ji -(<7i + 2a)). 



(B13) 
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5. Poles aX K = G + L — a 2 , K = G + L + <j 2 + 2a: These give two terms, 



7T cos 7ia 2 
sin7r(<T2 + a) 



-A 



2G+L-a 2 ,a 1 



G - a 2 , G + L - a 2 + a 

G + L + 1 - a 2 , G + 1 - 2 - a 



+ (<7 2 ->• -((7 2 + 2a)). 

(B14) 



These two terms are just the two terms in (B13) with o\ -H- o" 2 . 



Combining our results we have that 



n(L,U - \ -; ■ rA 2G+L _ a a2 1 

1 sm7r((Ti + a) 



Cr — 01, G + L — O"! + o; 

G + L + 1-0!, G + l-0i-a 



+ 3 sym 



+ cos(7ra) 



n=0 



(AG-Ti-Q,<TiAG+n+L+a,ff 2 + A G+n+ £ +Q , iCri A G _„_ Q>(T2 ) 

n + a, n + L + 2a 
n + 1, n + L + a + 1 



(B15) 



Here 3 sym refers to the three terms obtained by letting 01 — > — (01 + 2a), <7\ -H- 02, and 
0i ->■ -(02 + 2a), 2 ->■ 01. 

2. The G-sum 



Having computed the sum over K we have 



s ai a 2 {L) ■— ^ r 



G=0 



G + a, G + L + 2a 
G + l, G + L + a + 1 



(L; G) 



(B16) 



with H(L;G) given in (B15). First let us note that the infinite series in H(L;G) (see 
(B15)) gives a vanishing contribution when summed over G. The infinite series in H(L; G) 
contributes a term proportional to 



00 00 

EE 

G=0 n=0 



(Ac-n-a.o-i AG+ra+L+a,o- 2 + ^■G+n+L+a,a 1 h-G-n-a, 



<?1/ 



G + a, G + L + 2a, n + a, n + L + 2a 
G+l, G + L + a + 1, n + 1, n + L + a + 1 



(B17) 



Consider the change of variables G -H- n; under this action the gamma functions are invariant, 
as is A G+n+i+Q , )CT . . However, A G _„_ QjCr = — A„_ G „ Q , j(T , so in total the summand picks up a 
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(— 1) under the operation. As a result the double sum vanishes. This statement is true for 
all a e K. and a { e C. So we have that 



,{L) 



71 COs(7TO"i) 

sin7r(o"i + a) 



E A 



2G+L-cti,(T2 



G=0 



X r 



G + a, G + L + 2a, G - a u G + L - a ± + a 
G + 1, G + L + a + 1, G + L + 1 — &1, G + 1 — <7± — a 
+3 sym. 



(B18) 



We can now write S airT2 (L) in terms of four so-called "very well-poised" hypergeometric 
series (see A 3): 



(l) — ™ s l™iJ r 

sin 7r((Ti + a) 



Ar_ 



L — £7 1,(72 



X 7 V6 



a, —<Ji, L + 2a, L — o\ + a 
L + a + 1, L + 1 — <7i, 1 — o"i — a 

L — o"i — o" 2 L — <Ji + o" 2 + 2a 



L — o\ + a; a, — 01, L + 2a, 



+3 sym. 



(B19) 



These hypergeometric series have an excess of s = 4 — 4a, so are only absolutely convergent 
for a < 1. Assuming this, we may re- write the hypergeometric series using the transforma- 
tion (A13). The result is 



(L) 



7T COS(7RT]J 

2 sin n(ai + a) 



a, 2 — 2a, — a±, L + 2a, 2 + L — <j\ — a, 
1 — (Ti — a, L + a + 1, L + 1 — (Ti, 



L— (Ti — (72 L— (7i+(72+2a 
2 ' 2 

4+L— (7i— (72-4a 4+L-(7i+(72-2a 



1 + L — <j\ — a; 1 — a, 1 — a± — 2a, 1 + L, 



L — o\ — o" 2 L — (Ti + o" 2 + 2a 



+3 sym. 



(B20) 



3. Final result and checks 

In the previous section we computed the sum S ai<72 (L)] inserting this into (B2) yields 
(40). Let us examine the poles in p aia2 (L) as a function of a, <7j, and L. To do so it is 
useful to write this expression in terms of the regularized series 7V 6 (a] b, c, d, e, /) defined in 
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(A14): 



(7 1<T2 



xr 



1 COS^Ux) 

8n a sin7r((T 1 + a) 



2 — 2a, — <7i, L + 1, 1 + L — (Ti — a, 



L — <j\ — o" 2 L — o"i + o" 2 + 2a 



x 7 V 6 



1 + L — o\ — a; 1 — a, \ — U\ — 2a, 1 + L, 



2 2 

L — <j\ — a 2 L — <j\ + o" 2 + 2a 



+3 sym. 



(B21) 



The function iV§{a\ b, c, d, e, f) is entire in all its arguments, so the only possible poles arise 
from the gamma and trigonometric functions. 

In (B21) it appears that each of the four terms in have poles when a = 1/2, 1, 3/2, .... 
Upon inspection however one finds that the p CTl(T2 (L) is regular when a is a positive integer, 
and in these cases we may simplify our expression considerably. We record here the cases of 
a = 1,2: 

1 



XL) 



a=l 



16?r(L + 1) 



sin7r((7l + (72) -^f L - T (T2 )+(^yms) 



Sill 7TO"! Sin 7TO" 2 



+ 2n\, (B22) 



1 



a=2 



-1 
64^2 



L + l 
L + A 



x 



(L - <7i - <7 2 - 2)(L - (Ji + (x 2 + 2)(L + (7! - ct 2 + 2)(L + <n + o 2 + 6) 



x 



sin^ + ^ ^-a,-^ 

Sin 7TO"i Sin 7T(T2 



+ 2tt 



+ [(cot7T(7i)(<7i + 2)(L + 3) [L(L + 2) + <ti(<Ti + 4) - cr 2 ((7 2 + 4)] + {a x ^ a 2 )] |. 

(B23) 

For a = 1/2 the expression (B21) is also finite. However, for a — 3/2, 5/2, . . . the expression 
diverges. The divergences near a = 3/2 and a = 5/2 are given in(43) and (45). 

Next let us examine the pole structure of p ai(T2 (L) as a function of the mass parameters 
o"j. We restrict the ex, to be "on-shell", i.e. to have values corresponding to positive mass- 
squared (see section II). Under this restriction the only possible poles in (B21) are due to 
factor 1/ sin[7r((Tj + a)] and occur when <7j = —a + n, n G No- However, in the limit where 
<7j takes these values one finds that p Ul(72 (L) is regular. Thus, there are no poles in p Uia2 (L) 
as a function of (Jj. 
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Finally, let us examine the pole structure of p aiU2 (L) as a function of L. Recall that 
p aia2 (L) has been defined only for L G N . For these values it's clear that p aia2 (L) is 
regular. However, we may use (B21) to extend the definition of p ai(T2 (L) to L e C. In 
the complex L plane this expression has several possible poles arising from the Gamma 
functions 



L + 1, 1 + L — o\ — a, 



L — o"i — <7 2 L — o"x + o"2 + 2a 



(B24) 



2 ' 2 

and <7i permutations. For the poles at L = —1,-2,... one may explicitly compute the 
residues using Dougall's formula (A12); the residues of the four terms in (B21) cancel, so in 
fact p aiU2 (L) is regular for these values of L. Likewise, the Gamma functions r(l+L — — a) 
and permutations do not yield poles because their poles coincide with the zero of the series 
7Vq{cl] b, c, d, e, /) that occur when a is a negative integer. The remaining Gamma functions 
do indeed yield poles in p ai(T2 (L). We conclude that the expression (B21) has poles in the 
complex L plane at 

L = <j\ + er 2 — 2n, —<j\ + a 2 — 2a — 2n, <j\ — a 2 — 2a — 2n, —U\ — o" 2 — 4a — 2n. (B25) 

We may use Dougall's formula (A12) to compute the residue at these poles: 

r (ni -1 sin[7r(o-i + a 2 + a)] 

nes W J L=(T1+02 _2n 8 ir«T(a) sin tt^ + a) sinvr^ + a) 

1 — <7i — o"2 — 2a + 2n, n + a, n — ui, n — <r 2 , n — cti — o" 2 — a 



x r 



— <7i — o" 2 + 2n, 1 + 1 + n — a — <7i,l + n — a — <72, 1 



n 



o"i — <j 2 — 2a 



(B26) 



It is important to realize that (40) is not the unique extension of p Uia2 (L) to complex 
values of L. For example, under the assumption that L e N one may perform several 
hypergeometric transformations on (40) to derive an alternate expression for p aia2 (L) which 
agrees with (40) for L e No but has a different pole structure in the complex L plane. 



Appendix C: D = 3, 4 1-Loop corrections from 4-particle interactions 



Here we simply list the results of the calculations outlined in section IV A for diagrams 
shown in figure IV A for dimensions D = 3,4. The key point is that the constant A°" 2 (l) is 
given formally by setting Z = 1 in (24): 

l-d~ 



, „ n COS 7T (a + 

^ W 2 d+l 7T (d+3)/2 



-a, a + d, 



(CI) 
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This expression diverges for d — 1, 3, 5, . . . due to the factor T (1 — d/2). In these dimensions 
the divergence may be cancelled by the counterterm. 

For D = 3 the expression (CI) is finite. As a result, following the minimal subtrac- 
tion scheme (MS) we set the counterterms to zero and compute the self-energy correction 



ITi 



| D = g(l + a) cot na/(8ii£) , which represents a shift of particle l's mass M 2 — > Mf + IIi. 

For D = 4 — e we have the divergent expression 

(1 +a)(2 + a) 1 
e 



A'(l) 



D=A-e 



(1 + a)(2 + a) 



(C2) 



[-1 + 7 + % cot 7TO- - ln(47r) + 2^(3 + a)] + 0(e) 



16tt 2 £ 2 

where 7 is the Euler constant and ip{x) the digamma function. Defining the counterterms 

0(l + <7i)(2 + <7i) 1 



(5M1) 



g (I + a 2 )(2 + a 2 ) 1 
2 8tt 2 £ 2 e ' 



(5Mj) 



D=4-e 

leads to the self-energy correction 



D=4-e 



8n 2 £ 2 



(C3) 



n 



D=4-e 



32tt 2 



-1 + 7 + log 



47T 



+ 0(6) 



(C4) 



to the M 2 . As noted in section (IV A), these mass shifts encode the full context of the 
4-particle 1-loop corrections. Both of these expressions agree with the flat-space result 

9 r(i-f) 

flat 2(47r)^ 2 (M 2 ) 2 "^ • 



n 



(C5) 
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